A universal bound on iV-point correlations from inflation 
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Models of inflation in which non-Gaussianity is generated outside the horizon, such as curvaton 
models, generate distinctive higher-order correlation functions in the CMB and other cosmological 
observables. Testing for violation of the Suyama-Yamaguchi inequality tnl > (|/jvl) 2 , where /nl 
and tnl denote the amplitude of the three-point and four-point functions in certain limits, has been 
proposed as a way to distinguish qualitative classes of models. This inequality has been proved for a 
wide range of models, but only weaker versions have been proved in general. In this paper, we give 
a proof that the Suyama-Yamaguchi inequality is always satisfied. We discuss scenarios in which 
the inequality may appear to be violated in an experiment such as Planck, and how this apparent 
violation should be interpreted. We analyze a specific example, the "ungaussiton" model, in which 
leading-order scaling relations suggest that the Suyama-Yamaguchi inequality is eventually violated, 
and show that the inequality always holds. 
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A central goal of observational cosmology is to charac- 
terize the statistics of the initial perturbations of our uni- 
verse, thereby constraining the physics which generated 
these perturbations. Observations to date are consistent 
with adiabatic, Gaussian, scalar initial conditions, but 
deviation from a scale invariant initial power spectrum 
has been observed with 3er significance [1]. 

The search for non-Gaussian statistics in the initial 
fluctuations has emerged as a particularly interesting 
probe of inflation, due to the presence of many distinct 
signals (or "shapes" ) which probe inflationary physics in 
different ways. For example, detection of a nonzero three- 
point correlation function (£1^ £k 2 Ck 3 ) in "squeezed" con- 
figurations (fci <C min(A;2, A3)) would rule out all single- 
field models of inflation [2] , while detection of an "equi- 
lateral" three-point function (ki ~ fca ~ A3) is a generic 
test for self-interactions of the inflaton [3] . 

Throughout this paper, we use the following notation. 
We denote the connected part of an iV-point expecta- 
tion value by (-) c , and use the notation (•)' to denote 
the expectation value without the multiplicative factor 

(27r) 3 <5 3 (£ki). Fot a 3D field °"> we write A l( k ) = 
k 3 P a (k)/2TT 2 , where P„(k) is the power spectrum. 

The "local" model is a simple non-Gaussian model in 
which the initial adiabatic curvature is given by £ = 
Cg + §/Wl(Cg — (Cg))> where Cg is a Gaussian field and 
/jvl is a free parameter. The three-point and four-point 
functions in this model are: 
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type (1), but the relation between three-point and four- 
point functions is relaxed to an inequality: 



(CkiCk 2 Ck 3 }' — ■zfNL{Pk 1 Pk 2 
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eye.) (1) 
(CkiCksCkaCkJc = T Nh{PkiPk 3 P\\n 1+ -k 2 \ + 11 perm.) 



where P k = P ( (k) and t nl = (^Jnl) 2 - 

The local model can be generalized by introducing N 
Gaussian fields o\ , . . . , <jn (assumed for simplicity to be 
uncorrelated with equal power spectra) and taking the 
initial adiabatic curvature to be £ = Ai<n + Bjk(<Jj<7k — 
((TjCTfe)). This also gives rise to non-Gaussianity of local 
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tnl > \ ~^jnl 



(2) 



Suyama and Yamaguchi showed [4] that this inequality is 
always true at tree level in the (AN) expansion. In [5], it 
was observed that individual loop diagrams can violate 
the inequality, raising the interesting question of whether 
violation of the Suyama-Yamaguchi inequality could be 
an observational signature of loop diagrams in inflation. 

In this paper, wc will answer this question negatively: 
we will give a proof that the SY inequality is always sat- 
isfied. The proof is very general and does not depend 
on any physics; it comes from the requirement that a 
suitably constructed covariance matrix is always positive 
definite. We will also discuss scenarios in which the in- 
equality can appear to be violated observationally, and 
how this should be interpreted. 



GENERAL PROOF THAT t nl > (6f NL /5) 2 

The intuitive idea behind this proof is that if we define 
a local estimate of the small-scale power P, then tnl is 
the auto power of P on large scales, whereas (6/a/x/5) 2 
is the part of the auto power which can be attributed to 
the cross-correlation of P with £. 

To give the proof in maximum generality, we generalize 
the definitions of Jnl and tnl by taking the squeezed 
limit (k\ <C A2) of the three-point function and the col- 
lapsed limit (|ki + k2 1 <SC min(Aj)) of the four-point func- 
tion: 



/jVL = 



lim 



(CkiCk 2 Ck 3 



12 fe^o P c (fci)P c (A 2 ) 

(Ck 1 Ck 2 Ck 3 Ck 4 )c 



lim 



4 |k 1+k2 Ko P c (fci)P c (fc 3 )-Pc(l k i + k 2|) 



(3) 
(4) 
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This defines Jnl and tnl for arbitrary non-Gaussian ini- 
tial conditions. In the special case where the three-point 
and four-point functions have local shapes (1), we recover 
the usual definitions of Jnl and t nl . 

For wavenumbers k^ <C kg, let bg be a narrow band of 
wavenumbers near kg. Define a field Pk by: 



Pk 



1 

Vs 



d a k' Ck<k-k' 

k , e6s (27T)3 P(fc') 



(5) 



where Vg = Jy., eb d 3 k'/ (2tt) 3 is the volume of the band. 

The field P k represents the long- wavelength variation in 
the locally measured small-scale power. 

The power spectrum of Pk and its cross power spec- 
trum with (k are given by: 



(Pk Pk 



1 r d 3 kW (C- k < k '- ki Ck"C-k" + k i )' c 



v 2 
1 s 



-> 4rjvLPc(fcL) 



(2tt) 6 P c (k')P c (k") 
d 3 k' 2P c (fc')P c (fc") 
(2tt) 3 P c (fc')P C (fc") 
2 



(C^PkJ' = 4 



d 3 k' (C- k£ Ck< ki -k') 



(2tt)= 



P c (fc') 



12 



(6) 



(7) 



where "— >" denotes the leading behavior of each term in 
the Ul <C kg limit. 

Now consider the covariance matrix of the fields £, P: 



(Ck,p k *J' <4*A>' 



The determinant must be positive. 
Eqs. (6), (7), we get: 



TNL > 



5 fNL 



1 



2P c (k L )V s 



(8) 

Plugging in 
(9) 



Now take the limit (fci/fcs) — > 0. The second term on the 
RHS (which represents the disconnected part of the four- 
point function in Eq. (6)) goes to zero, and we obtain the 
SY inequality 



TNL > ( f /, 



TNL 



(10) 



This completes the proof. The proof is valid for ini- 
tial conditions with an arbitrary non-Gaussian PDF, and 
makes no assumptions about the physics. The SY in- 
equality emerges as a positivity constraint which is al- 
ways satisfied, in the same sense that the power spec- 
trum of any field is constrained to be nonnegative. The 
SY inequality was also interpreted recently as a positivity 
constraint in [6]. 



The proof above assumes that Jnl and tnl are de- 
fined in squeezed limits. If we define them at fixed scales 
k^,kg, then Eq. (9) shows that the SY inequality can be 
violated by an amount Atnl = l/(2Pf (fci)Vs). If wc 
assume that Vs = 0(k s ), then this means that sublead- 
ing contributions to the four-point function can violate 
the SY inequality, but such contributions must scale as 
Atjvl = 0(k\/k s ) in the collapsed limit. Given our 
definition of tnl, this corresponds to a four-point func- 
tion (CkiCk 2 Ck 3 Ck 4 ) which is finite in the collapsed limit. 
Conversely, a four-point function which is finite in the 
collapsed limit can generally have either sign and violate 
the SY inequality. Such examples can be interpreted as 
"accidental" tnl contributions from a four-point shape 
which is very different from the tnl shape. The amount 
by which the SY inequality can be violated is very small 
and difficult to detect with statistical significance in an 
experiment such as Planck, as we discuss next. 



ESTIMATORS 

We now study the question: if wc evaluate estimators 
for /nl and tnl in an experiment, do there exist realiza- 
tions of the data which appear to violate the inequality 
tnl > (I/wl) 2 , and if so, how should we interpret such 
a measurement? 

First consider an ideal experiment in which all modes 
£ k are measured without noise, over some range of scales 
in a periodic box. Given this data, we can write estima- 
tors for Jnl and tnl as follows. We normalize Fourier 
transforms so that (C k Ck') = Pc(k)VS^> , where V is the 
box volume. Define a field P k by: 



Pk = 



1 

Ns 



E 

k'e&s 



Ck'C: 



k-k' 



P ( (k') 



and estimators for Jnl and tnl by: 



Jnl 



tnl 



12V N L ^ P c (k) 



1 



4VN L 



keb L 



E 

kefci 



2V 2 



PkPk 



P C (fc) P c (k)Ng 



(11) 

(12) 
(13) 



Here, bL,bg are bins of wavenumbers near characteris- 
tic scales fci <C kg, and N a = X^kea 1 * s * ne number 
of modes in the bin (where a £ {L,S}). These defini- 
tions are closely analogous to the ones from the previous 
section, but here we are defining estimators which are 
applied to a single realization £ k , rather than taking an 
ensemble average over realizations. 

Since the sum of positive definite matrices is positive 
definite, we have the positivity constraint: 



DetE 



P c (k) 



CkCk 
CkPk 



C^Pk 

p k *pk 



> 



(14) 
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which gives: 



(15) 



where we have defined 
A 



V 



1/2 



2N L N S 



T — 

^ PA 



p c (k) 



(16) 
(17) 



Eq. (15) is the "estimator" version of the SY inequality. 
It applies to each realization £k individually (i.e. one docs 
not need to average over an ensemble of realizations for 
the inequality to apply). 

The quantity A which appears can be interpreted as 
a modification of the /ml estimator (/nl —> Inl/A) 
which is necessary for an estimator inequality to apply. 
Note that a similar modification was proposed in [7] as a 
way to reduce the estimator variance in the case where 
/nl is detected with statistical significance. 

The quantity Atnl is the maximum amount by which 
the estimator inequality tnl > (f/iVi/ A) 2 can be vi- 
olated by an individual realization. For the ideal ex- 
periment, Attvl is equal to the quantity l/{2P^{kL)Vs) 
obtained previously in Eq. (9). 

It is interesting to compare Atnl to the statistical 
error o~(tnl) which can be obtained using the estimator 
in Eq. (13). A short calculation shows: 1 



ct(t nl ) 



V 



E 



2^iVj^P c (fc) 2 



(18) 



1/2 

The ratio Atnl/o~ (tjv l ) is roughly N L ' . Assuming that 
the number of modes Nl is > 10, we conclude that in an 
ideal experiment, it is possible for a specific realization £k 
to violate the inequality tnl > (f/wi) 2 with statistical 
significance. 

How would we interpret such a violation if observed? 
As previously remarked, if we measure $nl and tnl at 
fixed scales fc/,, kg, the SY inequality can be violated by 
subleading terms which scale as tnl ~ 0(k L /kg). This 
scaling corresponds to a four-point function which is fi- 
nite in the collapsed limit. What we have shown here 
is that an ideal experiment can "see" such terms with 



1 This calculation makes the approximation that the connected 
four-point function of the field is zero on large scales. This 
type of approximation is common when forecasting higher-point 
estimators, e.g. in the context of CMB lens reconstruction, an 
analogous approximation has been shown to be accurate in [8]. 



statistical significance. In this scenario, the natural in- 
terpretation is that the tnl estimator is receiving acci- 
dental contributions from a four-point signal which is not 
the tml shape. 

The preceding discussion has considered an ideal ex- 
periment; let us now consider a more realistic case. For 
concreteness, consider a cosmic variance limited CMB 
experiment with £ max = 400. The statistical error 
&{ t nl)cmb for this experiment is « 4000. There will 
also be an estimator inequality of the form tnl > 
(§/jvlM) 2 - (Atjvl)cmb- Through numerical experi- 
ments (using a gradient minimization procedure to search 
for a realization which violates the SY inequality as 
much as possible), we find that (Atnl)cmb ~ 40000 
for this experiment. Since (Atjvl)cmb > (j(tjvl)cmb, 
it is possible to find CMB realizations in which the SY 
inequality appears to be violated with statistical signif- 
icance. For example, we can find realizations for which 
tnl = -40000 ± 4000, i.e. t nl is negative at 10u. Note 

1 /2 

that the ratio (Atnl)/o~{tnl) is roughly N L ' , where 
Nl is the number of "squeezed" CMB multipolcs with 
significant signal-to-noise; this number was found to be 
N L « 100 in [9]. 

When we apply a CMB estimator tnl to this dataset, 
we are measuring tnl averaged over a range of scales 
kL,ks, with a weighting which depends on CMB trans- 
fer functions (as opposed to our previously assumed 
flat weighting) but is peaked roughly at fcf, ~ 0.0004 
Mpc^ 1 and ks ~ 0.02 Mpc -1 . Now let us imag- 
ine we could do an ideal experiment in which we esti- 
mate tnl using all 3D modes £k throughout our Hub- 
ble volume, with the same weighting in kL,ks- In this 
ideal experiment, there will be a postivity constraint 

TNL > (I/jvlM) 2 - (ArjVL)ideal, where (ATjvx)ideal « 

(fci/fcsXATjvjJcMB = 800. 2 Since this constraint is sat- 
isfied for every realization £k, we do not actually need to 
do the ideal experiment in order to conclude that the in- 
equality tnl > —800 applies! Although there exist CMB 
realizations with tnl = —40000, it is mathematically im- 
possible to extend such a realization from the surface of 
last scattering to the Hubble volume; the extension au- 
tomatically satisfies tnl > —800. 

Given this picture, the only way to reconcile a mea- 
surement t^ b = -40000 ± 4000 (or more generally, a 
CMB measurement which violates tnl > (f/Wz,) 2 with 



2 To see that (Ar NL ) idcal rj {k L /k s )(^.T NL )cMB, we ar- 
gue as follows. For narrow bins in k^,kg, Eq. (17) gives 
(Ar]vx)idcal = (fc|_/2A 2 fc|)(A log kg). If we approximate the 
cosmic variance limited CMB as an ideal 2D measurement, 
then the value of Atnl will just be the 2D version of this, 
i.e. (Atjv_l)cmb = (fc 2 j /2A^fc|)(Alogfcs). This shows that 
(A-nvzJidcai « (k L /k s )(AT NL ) CMB for narrow bins; the gen- 
eral case follows by integrating over fc^ , kg with the appropriate 
weighting. 
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statistical significance) with positivity constraints seems 
to be to relax the assumption of translation invariancc. 
More precisely, translation invariance must be broken in 
a specific way where the four-point function near our 
surface of last scattering is very different from its mean 
value throughout the Hubble volume (so that tnl can be 
ps —40000 near the surface of last scattering, and have 
mean value > —800 in the Hubble volume). Since trans- 
lation invariance is assumed in the construction of tnl 
estimator, an observed violation of the SY inequality in 
the CMB would not be evidence for a particular model of 
inflation, but rather a sign that something is wrong with 
the assumptions which motivated searching for tnl in 
the first place. In this sense, testing the inequality with 
the CMB is more of a sanity check on the whole infla- 
tionary framework than a discriminator between models. 

The preceding analysis applies to any experiment 
where (At N l) > <t{t nl ) > (Ar^x) ideal- For exam- 
ple, this is the case for the upcoming Planck mission, 
since At nl /<7(t nl ) ps 10 and (At N l)/(At N l) ideal ps 
kh/ks ~ 300. In this case, it is possible to find CMB re- 
alizations in which the estimators violate the inequality 
tnl > (^Jnl/A) 2 with statistical significance, but it is 
mathematically inconsistent to interpret this as evidence 
for violation of the inequality tnl > (f Jnl) 2 throughout 
the Hubble volume. 



constant during inflation. This is a slight simplification 
of the scenario considered in [10] but contains all the 
qualitative features, including the t nl = O(f^f^) scaling 
at leading order in Jnl- 

The power spectrum in this model is: 



N 2 

P c (k)=N 2 P a (k) + ^P a2 (k) 



where 



P**(k) 



ri 3 q 

(2tt) 3 



P,(g)P CT (|k-q|) 



(21) 



(22) 



As a technical point, we note that if P a (k) is scale- 
invariant, then the power spectrum P a 2 (k) is infrared 
divergent. If the IR divergence is regulated by putting 
the fields in a finite box with length L, then the power 
spectrum diverges as P a 2(k) ~ 8tt 2 \og(kL)A^/k 3 . In 
what follows, we have left the regulator implicit, by writ- 
ing all IR-divergent quantities in terms of P a 2(k). 

A short calculation shows that the squeezed three- 
point function and collapsed four-point function are given 
by: 



(CkiCk 2 Ck 3 )' 
(CkiCk 2 Ck 3 Ck 4 )c 



N 3 

1 y (70 

2 

r4 



P CT 2(fc!)P CT (fc 2 ) 



(23) 



N^P a ,(\k 1 + k 2 \)P a (k 1 )P lT (k 3 ) 



UNGAUSSITON MODEL 

The so-called "ungaussiton" model from [10] is a mul- 
tifield model of inflation in which the three-point and 
four-point functions satisfy the following scaling relation 
in the limit of weak non-Gaussianity: 



r.A-2/3 ,4/3 
OZ V JNL 



(19) 



where C is a constant of order 1. Extrapolating to large 
/jvLj this scaling relation suggests that the inequality 
tnl > (|/wl) 2 is violated for f N L ^ A" 1 . 

As a check on our general theorem, in this section we 
will show explicitly that the Suyama-Yamaguchi inequal- 
ity is always satisfied in this model. We will find that the 
scaling relation (19) is only valid for Jnl *C A^ 1 , and 
takes a different form (which always satisfies the SY in- 
equality) for larger values of Jnl- Calculation of the 
three-point and four-point functions in the ungaussiton 
model also appeared recently in [11]. 

In the ungaussiton model, the initial adiabatic curva- 
ture is of the form 



l -N aa {a 2 (a 2 )) 



(20) 



where N^,N aa - are free parameters and the fields <f>, a 
are uncorrelated Gaussian fields with power spectra given 
by (fc 3 /27T 2 )P [T (fc) = (P//2tt) 2 , where Hj is the Hubble 



Plugging the above expressions into the defini- 
tions (3), (4) of /atl and tnl, we find: 



6 a 3 /3(k L ) 

5 JNL (l + a 2 /?(fc L ))(l + a 2 /3(fc s )) 



tnl 



a 4 (3(k L ) 



(l + « 2 /3(fc L ))(l + a 2 /3(fc s )) 2 



(24) 
(25) 



where we have defined dimensionless parameters a = 
N aa /N 2 and 0(h) - N 2 P a 2 (k) / (AP a {k)) . 

Note that we recover the scaling relation (19) in the 
Gaussian limit (i.e. a 2 (3 <C 1 with f3 = C(A 2 )). However, 
in the general form in Eqs. (24), (25), it may be seen that 
the SY inequality is always satisfied. Indeed, 



tnl 



(6W5) 2 



1 



1 



u 2 fi{k L 



(26) 



and the RHS is always > 1 since p{k£) must be positive. 
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